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Thermotropic interface and core relaxation dynamics of 
liquid crystals in silica glass nanochannels: A dielectric spectroscopy 
study 


Silvia Cahis," Lech Borowik," Andriy V. Kityk,*“ Manfred Eich, ^ Mark Busch‘S and Patrick Huber* 


We report dielectric relaxation spectroscopy experiments on two rod-like liquid crystals of the cyanobiphenyl family (5CB 
and 6CB) confined in tubular nanochannels with 7 nm radius and 340 micrometer length in a monolithic, mesoporous silica 
membrane. The measurements were performed on composites for two distinct regimes of fractional filling: monolayer coverage 
at the pore walls and complete filling of the pores. For the layer coverage a slow surface relaxation dominates the dielectric 
properties. For the entirely filled channels the dielectric spectra are governed by two thermally-activated relaxation processes 
with considerably different relaxation rates: A slow relaxation in the interface layer next to the channel walls and a fast relaxation 
in the core region of the channel filling. The strengths and characteristic frequencies of both relaxation processes have been 
extracted and analysed as a function of temperature. Whereas the temperature dependence of the static capacitance reflects 
the effective (average) molecular ordering over the pore volume and is well described within a Landau-de Gennes theory, the 
extracted relaxation strengths of the slow and fast relaxation processes provide an access to distinct local molecular ordering 
mechanisms. The order parameter in the core region exhibits a bulk-like behaviour with a strong increase in the nematic ordering 


just below the paranematic-to-nematic transition temperature 
surface ordering evolves continuously with a kink near Tpf^. 
with the complete filling reveals that the molecular order in tht 
molecular layers at the wall. 

1 Introduction 

Composites of liquid crystals (LCs) and monolithic, meso¬ 
porous solids, optically transparent porous silica in particu¬ 
lar, are promising hybrid materials for organic electronics and 
the emerging field of nano photonics^"-. They can be eas¬ 
ily prepared by melt infiltration, profit from the mechanical 
stability of the porous solid template and from the large varia¬ 
tion in electrical and optical properties offered by the plethora 
of different liquid crystalline/mesoporous host combinations 
nowadays available. Moreover, they allow one to explore the 
thermodynamics, structure and transport characteristics of liq¬ 
uid crystalline systems in restricted geometries, and thus phe- 
nomenologies which are of high interest both in nanoscience 
and nanotechnology-. 

Most prominently, the paranematic phase has been widely 
explored in theoretical-"— and experimentali^— studies at 
planar surfaces and in porous media. It is characterised by a 
residual nematic order and thus the absence of a ’’true” isotropic 
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TpN and subsequent saturation during cooling. By contrast, the 
A comparison of the thermotropic behaviour of the monolayer 
core region of the pore filling affects the order of the peripheral 


liquid state at high temperatures. The evolution of the orien¬ 
tational order parameter from this pre-ordered state to the ne¬ 
matic phase can be described by a ’’nematic ordering field”, a 
within a Landau-de-Gennes free energy approach-"—. The 
strong first order I-N transition is replaced by a weak first 
order or continuous paranematic-to-nematic (P-N) transition 
at a temperature Tp^ and may also be accompanied by pre- 
transitional phenomena in the molecular orientational distri¬ 
bution—. For tubular pore geometry, this effective field is 
strongly dependent both on the average pore radius R(oc R~^) 
and on the LC-wall interaction. 

Optical polarimetr y 17^20^22^23 pj-ovides arguably the most 
accurate insights in the order parameter behaviour in the vicin¬ 
ity of the paranematic-to-nematic transitioniSiii. Note, how¬ 
ever, that this technique and most other experimental methods 
probe an effective (averaged) molecular ordering. 

By contrast computer simulations on LCs in thin film and 
pore geometry can give spatially resolved informatio n 1242 ^ 13 ^ 24 - ^ 
These studies indicate pronounced spatial heterogeneities, en¬ 
compassing interface-induced molecular layering and radial 
gradients both in the orientational order and reorientational 
dynamic s 

Dielectric relaxation spectroscopy has proven as particu¬ 
larly suitable to achieve such local information on confined 
LCsi^"— . A number of studies on LCs in porous media evi- 
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Fig. 1 Sketch of the capacitor geometry for the dielectric relaxation 
experiments: Gold electrodes were evaporated onto the mesoporous 
silica membrane, so that the sample forms a simple parallel circuit 
both in the regime of monomolecular layer coverage (a) and for the 
completely filled channels (b). (c) Illustration of two regions with 
two distinct molecular mobility for completely filled nanochannels: 
the interface layer next to the channel wall with a slow dipolar 
relaxation and the core region characterised by fast relaxation 
dynamics. 


dently document that the rate of dipolar relaxations (and thus 
orientational and translational mobility) usually differ signif¬ 
icantly between the molecules in the pore wall proximity and 
the ones in the channel centre^"—. 

In dielectric studies Aliev et al.— — found a slow surface 
mobility in comparison to the dynamics in the pore centre for 
liquid crystals confined in tortuous and tubular mesopores. 
Moreover, a significant broadening of the dielectric spectra 
was observed and traced by the authors to inhomogeneous 
couplings of the molecules to the pore walls and coupling vari¬ 
ations among the molecules themselves^. 

It turned out that the rate of dipolar relaxation differs by 
about two orders of magnitude between the molecules located 
at the host-guest interface and the ones located in the core re¬ 
gion of the pore filling, see sketch in Fig. Ic. Therefore, both 
relaxation processes can be discriminated and the correspond¬ 
ing dipolar relaxation strengths provide then access to local 
molecular ordering. Employing this distinct dynamical be¬ 
haviour, we could recently show for 7CB, a prominent mem¬ 
ber of the rod-like cyanobiphenyl nematogens that the order¬ 
ing in the core region is reminiscent of the bulk behaviour—. 

It is characterised by an abrupt increase in the nematic order¬ 
ing in the transition region. By contrast, the surface ordering 
exhibits a continuous thermotropic evolution of nematic order 
with a gradual change in slope and an asymptotic vanishing in 
the paranematic phase. 

In this article we extend our previous dielectric spectroscopy 
studies on two other members of the cyanobiphenyl family, 
i.e. ^CB (^=5,6) embedded in silica membranes with par¬ 
allel aligned nanochannels. Overall, we find a very similar 
behaviour of the effective average as well as local order pa¬ 
rameters as in the case of 7CB. As we will outline below, it is 
again possible to describe the experimental observations in a 
phenomenological manner by applying a Landau-De Gennes 
model. The analogous findings and semi-quantitative descrip¬ 
tions document that the dielectric spectroscopy technique along 
with the phenomenological approach and key physical princi¬ 
ples outlined in Ref.— seem to be applicable to confined ne¬ 
matics in general. 

2 Experimental 

The nematic LCs 5CB and 6CB have been purchased from 
Merck, Germany. The porous silica (/7Si02) membranes were 
obtained by electrochemical anodic etching of highly /7-doped 
(100) silicon wafers which have been subjected to thermal ox¬ 
idation for 12 h at r=800 under standard atmosphere. The 
resulting array of channels are aligned along the [100] crystal¬ 
lographic direction, i.e. perpendicular to the membrane sur¬ 
face. The average channel radius is R = 6.6±0.5 nm (poros¬ 
ity P = 55±2%) as determined by recording volumetric N 2 - 
sorption isotherms aiT=11 K. For the dielectric measurements. 
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gold electrodes have been deposited onto the porous mem¬ 
brane. All measurements were performed on samples cut from 
a monolithic porous membrane of d =340 jim thickness. Two 
samples with electrode area of 97 mm^ (geometric capacitance 
Co = EoS/d = 2.52 pF) and 104 mm^ (Co = 2.71 pF) have been 
filled by nematic LCs 5CB and 6CB, respectively. 

The entirely filled samples (fractional filling / =1.0) have 
been obtained by capillary imbibition of the LC melt—. To 
prepare the partially filled samples with layer coverage we im¬ 
mersed it into a binary LC/cyclohexane solution (5 vol. %) for 
about 20-30 minutes. After evaporation of the high-vapour- 
pressure solvent (cyclohexane), the low-vapour-pressure LC 
remained in the porous matrix. The final filling fraction / and 
the complete evaporation of the solvent have been verified by 
comparing the sample weight before and during the adsorption 
procedure^. The obtained fraction fillings / = 0.16 ±0.01 
(for 5CB in pSi02) and / = 0.18 ± 0.01 (for 6CB in pSi02) 
correspond approximately to one monomolecular layer. 

Dielectric spectra have been recorded in the frequency range 
from 0.5 kHz to 15 MHz using the impedance/gain-phase anal¬ 
yser Solatron-1260A. The measurements have been performed 
at selected temperatures between 292 K and 324 K and con¬ 
trolled with an accuracy of 0.01 K by a LakeShore 340 tem¬ 
perature controller. 

Note that both, the entirely filled samples and the partially 
filled samples with monolayer coverage form simple parallel 
electrical circuits, as schematically sketched in Figs, la and 
lb, respectively. The advantage of the simple geometry is ev¬ 
ident: Since the channels are all arranged parallel to the ex¬ 
ternal electric field, the effective complex permittivity of the 
composite, £*, is given by the permittivities of the compo¬ 
nents ^si 02 (silica substrate), (nematic LC) and Cy = 

1 (vacuum permittivity) weighted with the corresponding vol¬ 
ume fractions: e* = (1 ±F[/£*^5±(1 —f)£y]. In the 

frequency range 0.0005-15 MHz £^-q^ ^ £si 02 
not exhibit any frequency dispersion. Thus, the relaxation 
behaviour observed in the effective permittivity, e*(v) corre¬ 
spond to that of the confined LC. For more complicated pore 
geometries, e.g. for partial fillings in the capillary condensa¬ 
tion regime depolarisation effects may lead to phase shifts of 
the internal electric field compared to the external field, which 
considerably complicates the analysis - see also Ref.—. Of 
course, the channel walls of the present matrix are rough, and 
the radius is likely to vary along the long channel axis which 
may result in deviations from the ideal parallel geometry. Nev¬ 
ertheless, the parallel-circuit assumption is expected to be a 
good approximation. 


3 Results and discussion 

3.1 Molecular mobility probed by dielectric spectroscopy 

In Fig. 2 we display the frequency dispersion of the imaginary 
capacitance C"(v) = e"(v)Co (symbols) at five selected tem¬ 
peratures, r, for entire and partially filled nanoporous silica 
membranes. 

In order to analyse the observed relaxation behaviour in 
more detail, we resort to a representation of the dielectric data 
in so-called Cole-Cole diagrams, i.e. we plot the imaginary 
part C"(v) on the vertical axis and the real part C"(v) on the 
horizontal axis with frequency V as the independent param¬ 
eter, see insets in Fig. 2. In such a diagram, a material that 
has a single relaxation frequency, as typical of the classical 
Debye relaxator, will appear as a semicircle with its center ly¬ 
ing on the horizontal at C = 0 and the peak of the loss factor 
occurring at 1 /t, where T is a measure of the mobility of the 
molecules (dipoles). It characterises the time required for a 
displaced system aligned in an electric field to return to 1 of 
its random equilibrium value (or the time required for dipoles 
to become oriented in an electric field). A material with mul¬ 
tiple relaxation frequencies will be a semicircle (symmetric 
distribution) or an arc (nonsymmetrical distribution) with its 
center lying below the horizontal at C = 0. In that sense, the 
Cole-Cole representation allows one to geometrically illus¬ 
trate and analyse the relaxation behaviour of a given system 
in a quite simple manner.^i^ An analysis of the experimen¬ 
tal data shows that the dielectric dispersion of the complex 
capacitance C*{co) can be well described by two Cole-Cole 
processes: 

C*((o) = e*(m)Co = 

ACi AC2 

^^l±(imTi)i-«i l±(imT2)i-«2’ 

where m = 2;rv is the cyclic frequency. Coo = £ooCo is the 
high frequency limit capacitance expressed via the high fre¬ 
quency permittivity Coo, ACi = AEiCq and AC 2 = Ae 2 Co are 
the capacitance relaxation strengths expressed via the dielec¬ 
tric relaxation strengths Aei and Ae 2 of the slow process I and 
the fast process II, respectively, Ti and T 2 are the mean relax¬ 
ation times of the corresponding processes. The Cole-Cole 
exponents ai and a 2 (0 < ai,a 2 < 1) characterise the dis¬ 
tribution of the relaxation times; its zero limit corresponds 
to the Debye process with a single relaxation time. Solid 
lines in Fig. 2 are the best fits as obtained by a simultane¬ 
ous analysis of the measured real and imaginary capacitance 
based on Eq. 1. The deviations of the experimental data points 
from the fitting curves at low frequencies originate in ionic dc- 
conductivity. For this reason the low-frequency region was al¬ 
ways excluded from the fitting analysis.This fitting procedure 
has been applied to the measured dispersion curves in the en- 
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Fig. 2 Frequency dispersion of the imaginary capacitance C"(v) in 5CB-/7Si02 and 6CB-pSi02 nanocomposites for five selected temperatures. Panels (a) and 
(c) correspond to partially filled matrices 5CB-j!7Si02 (/ = 0.16) and 6CB-/7Si02 (/ = 0.18), respectively, in the regime of monomolecular layer coverage. 
Panels (b) and (d) correspond to nanocomposites with entirely filled pores (/ = 1.0). Insets in panels (a)-(d) displays Cole-Cole plots for three selected 
temperatures. Symbols in the panels and inserts are the measured data points. Solid lines are the best fits based on the Cole-Cole relaxation model (see Eq.(l)). 
The geometric capacitance, Cq, equals 2.52 pF (5CB-pSi02) and 2.71 pF (6CB-pSi02). 
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tire temperature range. In the case of entirely filled samples, 
5CB-/7Si02 and 6CB-/7Si02 (/ = 1.0), uncertainties regard¬ 
ing the extracted fit parameters can occur at higher tempera¬ 
tures, particularly above the paranematic-to-nematic transition 
point, Tp^. Then the maximum of the imaginary part, C"{v), 
shifts out of the upper limit of the frequency window (v^ > 15 
MHz) employed in the dielectric measurements. Thus only a 
part of the left wing of the relaxation band is observed. 

Fortunately, this problem can be resolved by following the 
route described in detail in Ref.—, where one takes into ac¬ 
count that the capacitance relaxation strength AC 2 can be de¬ 
termined alternatively by employing Eq. (1) in its static limit 
(V ^ 0): 

AC 2 = Cst — Coo — ACi. (2) 

Here the static dielectric constant, Cst, is determined di¬ 
rectly from the Cole-Cole plot by its extrapolation to the low 
frequency region. Obviously, Cst is practically independent 
from other extracted fit parameters. Coo, on the other hand, 
it is strongly dominated by the electronic polarisability which 
is weakly temperature dependent. Estimations show that its 
contribution to the temperature changes of the static capaci¬ 
tance, similarly as in the case of 7CB-/7Si02— , do not exceed 
5% in the entire T-range, i.e. it remains within the error of 
the fitting analysis. In the fitting procedure the best parameter 
sets have been determined by a minimalisation of the differ¬ 
ence between the AC 2 values extracted by the two alternative 
ways: (i) direct fits of measured dispersion C'(v) and C"(v) 
and (ii) extracting Cst-value from the Cole-Cole plots and sub¬ 
sequent employing of Eq.(2). The combination of these pro¬ 
cedures ensures unambiguity of the extracted fit parameters. 
The corresponding temperature dependences are displayed in 
Figs. 3-6. 

The interpretation of the dispersion curves is similar to that 
reported in Ref.— as well as in a series of previous dielec¬ 
tric studie s 39^41-43^45^56^57 ^ entirely filled samples (/ = 

1) the slow relaxation corresponds to the rotational dynam¬ 
ics of the molecules in direct contact with the pore walls (so- 
called surface or more strictly spoken interface relaxation), 
see sketch in Fig. Ic. The fast relaxation process originates 
in the molecular rotational dynamics in the core region (the 
so-called core relaxation). It is obvious that for partially filled 
samples, i.e. in the regime of monomolecular layer coverage, 
the dielectric relaxation is strongly dominated by the surface 
relaxation. 

Nevertheless, a weak contribution of the fast relaxation is 
also present. Its presence is indicated by the slightly asymmet¬ 
ric shape of Cole-Cole plots observed at higher temperatures. 
The corresponding relaxations strength, AC 2 (see Fig. 6, / = 
0.16 (5CB) and / = 0.18 (6CB)), decrease with decreasing 
temperature and practically vanish below 290 K. We believe 
that due to increasing intermolecular spacings at higher tem¬ 
peratures some molecules are pushed into the next molecular 



T[K] 


Fig. 3 Cole-Cole parameters of the slow relaxation process I, a\ , and fast 
relaxation process II, 0^2 vs T as extracted within the fitting procedure for 
5CB-/7Si02 (a) and 6CB-/7Si02 (b) nanocomposites in the regimes with 
entirely (/ = 1.0) and partial (/ = 0.16 (5CB-/7Si02); / = 0.18 
(6CB-/7Si02)) filled channels. 

layer, resulting in a considerably faster dipolar relaxation. The 
’’fast” relaxation process in the monomolecular layer regime, 
on the other hand, differs qualitatively from the one observed 
in the core region of entirely filled matrices. In the first case 
the molecules are located at the tbp interface. In the sec¬ 
ond one they are located in the core region of the pore filling 
and, accordingly, are infiuenced by the collective molecular 
ordering resulting from the paranematic-to-nematic transition. 
Therefore, the smooth temperature variations of all Cole-Cole 
parameters observed in the monomolecular layer regime con¬ 
trast with the ones observed for the entirely filled matrices. 
Roughly speaking, the monomolecular layer regime does not 
show a phase transition. This is evident if one compares e.g. 
the temperature dependences of the static capacitance, Cst(T), 
for both regimes of pore filling, see Fig. 5. 

For the nanocomposites with entirely filled matrices (/ = 1.0) 
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Fig. 4 Relaxation times of the slow relaxation process I, Ti, and fast 
relaxation process II, T 2 vs 7“^ as extracted within the htting procedure for 
5 CB-/ 7 Si 02 (a) and 6 CB-pSi 02 (b) nanocomposites in the regimes with 
entirely (/ = 1 . 0 ) and partial (/ = 0.16 ( 5 CB-pSi 02 ); / = 0.18 
( 6 CB-/ 7 Si 02 )) filled channels. 


the Cole-Cole parameters, ai and a 2 (see Fig. 3) exhibit op¬ 
posite tendencies in their temperature evolution. The surface 
relaxation (process I) is characterised by a rather broad distri¬ 
bution of the relaxation times at room temperature (ai ^ 0.48- 
0.50), but becomes somewhat narrower at high temperatures 
and approaches a magnitude of 0.35-0.40 in the paranematic 
phase. The dipolar relaxation in the core region (process II), 
on the other hand, is characterised by a narrow distribution of 
the relaxation rates at room temperature (a 2 ^ 0.05-0.06) and 
rises to about 0.12-0.14 in the paranematic state, with a char¬ 
acteristic, step-like change in the vicinity of the paranematic- 
to-nematic transition. Similar changes in the phase transition 
region are observed in the temperature evolution of other re¬ 
laxation parameters, T 2 (T) (Fig. 4) and AC 2 {T) (Fig. 6). This 
means that an orientational ordering in the core region affects 
the relaxation parameters. For comparison, the surface relax¬ 
ation parameters, ai, Ti and ACi exhibit here only smooth 
temperature variations. At this stage, one can already con¬ 
clude that the molecular ordering in the core and surface re¬ 
gions of the pore filling are expected to be considerably dif¬ 
ferent. This issue will be discussed in details below. 

For the pores with monomolecular layer coverage, the pa¬ 
rameter ai exhibits quite similar temperature variation as the 
one found for the entirely filled samples. The corresponding 
ai (r)-curves are slightly shifted down, see Fig. 3. 

The relaxation times of the surface and core processes are 
displayed in Fig. 4. The fast relaxation in the core region ex¬ 
hibits Arrhenius-like behaviour (T 2 = Toi^wi^ai/ktT)) with 
somewhat different activation energies in the nematic (T « 
Tpn) and paranematic (T » Tp^) phases. An exception is 
the vicinity of the paranematic-to-nematic transition, where 
smeared, step-like variations are observed. They may be at¬ 
tributable to a specific behaviour of the attempt relaxation time, 
To2^ presumably because of changes in the activation entropy 
caused by the orientational molecular ordering at the phase 
transformation. The surface relaxation, Ti, on the other hand, 
exhibits only a smooth temperature variation accompanied by 
a gradual change in slope, but with no special features in the 
vicinity of Tp^. This process is much slower because of the 
large viscosity in the surface layer compared to the bulk vis¬ 
cosity. This hinders the rotational dynamics of the moleculesi^. 
Both relaxation times rise with increasing orientational order¬ 
ing which is typical of rotational dynamics around the short 
molecular axis^^iS. 


3.2 Thermotropic nematic order: Landau-de Gennes Anal¬ 
ysis and Surface/Core Partitioning 

In the Landau-de Gennes theory the nematic (orientational) or¬ 
der parameter is defined as 2 = ^ (3 cos^ O — l), where 0 is the 
angle between the axis of the molecule and the director, H, and 
brackets mean an averaging over the ensemble of molecules. 
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Fig. 5 Static capacitance, Cst as extracted within the fitting procedure 
(symbols) for 5CB-pSi02 (/ = 1.0(a), / = 0.16(b)) and 6CB-/7Si02 
(/ = 1.0(c), / = 0.18(d)) nanocomposites. Solid blue lines are the best fits of 
Cst (r)-dependence (/ =1.0) based on the KKLZ approach. Dash-dot blue 
lines indicate isotropic baselines, (T), obtained within the same fitting 
procedure. 


The size of an ensemble domain is an important issue in such a 
description, particularly when the order parameter is spatially 
inhomogeneous. Smaller sizes are appropriate to describe lo¬ 
cal properties of the molecular ordering. Larger ones are more 
suitable in a characterisation of an effective (averaged) value 
of the order parameter. Its appropriate choice is thus a matter 
of the specific system under consideration and/or the approach 
used for its description. Specifically, the surface layer and core 
region of the systems studied here may be considered as two 
regions with homogeneous order behaviour, characterised by 
two distinct, local order parameters, Qs and Qc, respectively. 
The effective (averaged) order parameter over the pore filling, 
2, can then be expressed as a superposition of the elementary 
contributions: Q = where Qs = WsQs and Qc = WcQc, 

and Wc (w^ + = 1) are the weight factors (volume frac¬ 

tions) of surface and core components, respectively. Although 
quite simplistic, this approach has an evident practical bene¬ 
fit: All three quantities, i.e. Q, Qs and Qc can be extracted 
independently from the experiment. Following the approach 
developed in Ref./^i, which is based on the Maier and Meier 
equation— the excess changes of the static capacitance are 
proportional to Q: 

Cst-4""^oc2, (3) 

where is the isotropic baseline of the static capaci¬ 
tance. In the case of the bulk nematic LC (T) represents 

the bare temperature dependence of the capacitance relaxation 
strength in the isotropic phase, whereas its extrapolation to the 
nematic phase gives an isotropic baseline relative to which an 
excess contribution due to an orientational ordering (/7-term) 
is counted. Similarly, an excess changes of the capacitance re¬ 
laxation strengths ACi and AC 2 are proportional to Qs and Qc, 
respectively: 

A.Ci{T) - o. Q^{T)- (4) 

AC2(7’) - Acf‘'^(7’) oc Qc{T). 

where {T) and AC 2 (r) are the corresponding isotropic 

baselines. The factor of proportionality depends on the molec¬ 
ular dipole moment and its orientation with respect to the long 
principal axis of the molecule, the molecular number density, 
the internal field factors and the temperature, but it is identical 
in Eqs. 3 and 4. 

In Fig. 5 we display the temperature dependences of the 
static capacitance for the composites 5CB-/7Si02 and 6CB- 
/7Si02. In the monomolecular layer regime [5CB-/7Si02 (/ = 
0.16), section (b); 6CB-/7Si02 (/ = 0.18) section (d)] only a 
gradual increase of Cst is observed upon cooling indicating a 
weak orientational ordering, but no hints of a phase transfor¬ 
mation. For the entirely filed sample [5CB-/7Si02 (/ = 1.0), 
section (a); 6CB-/7Si02 (/ = 1.0) section (c)], on the other 


This journal is ©The Royal Society of Chemistry [year] 


Journal Name, 2010, [vol], 1 - TT^ | 7 
















Fig. 6 Capacitance relaxation strengths of the slow relaxation process I, 
ACi and fast relaxation process II, AC 2 vs T as extracted within the fitting 
procedure for 5CB-pSi02 (panel (a), / = 1.0, / = 0.16) and 6CB-/7Si02 
(panel (b), / = 1.0 / = 0.18) nanocomposites. Dash-dot blue lines indicate 
the baselines of the isotropic state, see labels. 


hand, pronounced changes with a characteristic kink at Tpj^ on 
C^i{T) dependences are observed for both nanocomposites. A 
typical feature of these dependences is their continuous char¬ 
acter with a precursor behaviour quite similar to that observed 
in recent optical birefringence studiesi2i22i2Si2^. 

Whereas the molecular ordering in the bulk, i.e. the isotropic- 
to-nematic phase transition, can be described by a Landau- 
de Gennes theory, Kutnjak, Kralj, Lahajnar, and ZumeriSaii 
have extended the corresponding phenomenological approach 
(hereafter denoted as KKLZ model) towards nematic phases 
in cylindrical confinement. In principal, this model is appli¬ 
cable to entirely filled channels, only and its description de¬ 
pends on the anchoring conditions at the channel walls. The 
untreated silica channels enforce planar anchoring— with a 
preferred orientation of LC molecules along the long channel 
axes. A key point of the KKLZ theory is a bilinear coupling 
between the order parameter and the nematic ordering field, 
cr, which results in a residual nematic ordering (paranematic 
state) for T > instead of an isotropic one. Omitting a 
description of the KKLZ model and details of the fitting pro¬ 
cedure, which can be find in Refs.Ji22a2Sa22, we present here 
only the main results of this analysis. Figs. 5a and 5c dis¬ 
play the best fits of measured Cst(T)-dependences along with 
the isotropic baselines, (T), obtained within the same fit¬ 
ting procedure. The difference between the bulk and effective 
transition temperatures, AT*, which calibrates the temperature 
scale of the KKLZ-model, was taken to be equal 3.4 K (5CB) 
and 3.2 K (6CB) based on the results of recent optical polari- 
metric studies^. The best fits yield cr-values of 0.91 (5CB- 
pSi 02 ) and 1.10 (6CB-/7Si02), which corresponds to a super¬ 
critical regime in both cases. The magnitudes of the critical 
radius, Rc = 2Rg equal to 12.0 nm (5CB-pSi02) and 14.5 nm 
(6CB-/7Si02), respectively, which within an experimental er¬ 
ror (± 1.0 nm) is in agreement with that obtained in the optical 
birefringence studies: 12.1 nm and 14.0 nm— , respectively. 
Obviously, the KKLZ approach gives an adequate description 
of the paranematic-to-nematic transition. Deviations of the 
experimental data points from the fit curves observed at lower 
temperatures originate in order parameter saturation, which is 
not appropriately described by the free energy expansion of 
the KKLZ-model, since it is limited to a fourth-order expan¬ 
sion of the free energy. 

In accordance with Eq.(3) the difference, Cst(T) — {T) 

is proportional to Q{T). Normalised to a Q-value at T = 
290 K, the 2(T)-dependences of 5CB-pSi02 and 6CB-;?Si02 
nanocomposites are displayed in the upper panels of Figs. 7a 
and 7b, respectively. Using the ACi(r) and AC 2 {T) depen¬ 
dences presented in Fig. 6 for entirely filled matrices (/ = 1.0) 
and Eqs.(4) the effective (averaged) order parameter Q can 
be deconvoluted into elementary contributions characterising 
molecular orderings in the surface (Qs) and core (Qc) regions. 
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Fig. 7 The effective (averaged) orientational order parameter, QysT and its deconvolution in the elementary contributions (averaged local order parameters) 
describing molecular orderings in the surface (Qs(T)) and core (Qc{T)) regions for 5CB-pSi02 (/ = 1.0, section (a)) and 6CB-/7Si02 (/ = 1.0, section (b)) 
nanocomposites. For convenience all order parameters are normalised to the value Q at 290 K. The relation Q — Qs-\-& holds for each T. 
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A principal challenge are the determination of the isotropic 
baselines, Acf (T) (/ = 1,2). We trace first Acf (T). Above 
TpN AC 2 {T) rather fast saturates to a nearly temperature inde¬ 
pendent value. A linear extrapolation of this dependence be¬ 
low TpN, as it is displayed in Fig. 6 represents the isotropic 
baseline in the confined nematic phase. This approach is not 
applicable, however, to describe the (T)-baseline. ACi (T) 

exhibits evident nonlinear asymptotic behaviour extending far 
above Tp^. Fortunately, having two other isotropic baselines it 
can be calculated as: Acf“^(r) =cf“^(r)-Coo-Acf°^(r), 
see labeled dash-dotted lines in Fig. 6 . The difference AQ (T) — 
A^(iso)(y.) _ ^ 2 ^ jg proportional to elementary contribu¬ 

tions, Qs and Qc. Normalised again to Q-value at (T =290 
K), 2^(7) and Qc{T) dependences are displayed in the lower 
panels of Fig. 7a and 7b as for 5 CB-pSi 02 (/ = 1.0) and 6 CB- 
/ 7 Si 02 (/ = 1.0), respectively. Variations Qs{T) and Qc{T) in 
the region of Tp^ are considerably different. 

The molecular ordering in the core region is reminiscent of 
the evolution of the bulk order parameter. It exhibits a strong 
increase in the nematic ordering just below Tp^ and its sub¬ 
sequent saturation at T << Tp^. The surface ordering, on the 
other hand, exhibits a continuous change with a smeared kink 
near Tp^ and an asymptotic decrease above this temperature. 
Moreover, the molecular ordering in the core region has an 
influence on the ordering in the periphery, particularly in the 
molecular layer located next to pore walls. This can be in¬ 
ferred by comparing the ACi(r)-dependences measured for 
the complete and partial pore fillings, see Fig. 6 . Particularly, 
a change in slope below Tp^ is quite obvious only for the sam¬ 
ples with entirely filled pores. Presumably, the molecular or¬ 
dering in the core region is transferred to the periphery via 
intermolecular interactions. 

4 Conclusion 

In conclusion, we reported a dielectric study on the calamatic 
nematic crystal 5CB and 6 CB confined in cylindrical nanochan¬ 
nels of monolithic silica membranes. The measurements have 
been performed on composites with two distinct regimes of 
pore fillings: monomolecular layer coverages and entirely filled 
pores. Whereas for the composites with monomolecular layer 
coverages a slow surface relaxation dominates the dielectric 
properties, the dielectric spectra of nanocomposites with en¬ 
tirely filled pores can be well described by two Cole-Cole 
processes with well separated relaxation times. The fast re¬ 
laxation is associated with a rotational dynamics of molecules 
in the core region of the pore filling. The slow relaxation orig¬ 
inates from a surface LC layer next to the pore walls. 

In the regime of monomolecular layer coverage the tem¬ 
perature evolution of the static capacitance exhibits a smooth 
temperature behaviour with no hints of a phase transformation 


in the entire temperature region. This is in contrast to the ma¬ 
trices with entirely filled pores. The static capacitance of such 
samples evidently demonstrates anomalous behaviour with a 
precursor behaviour caused by orientational molecular order¬ 
ing due to the paranematic-to-nematic transition. For the en¬ 
tire pore filling it can be well described by a phenomenological 
Landau-de Gennes theory (KKLZ approach). The correspond¬ 
ing analysis yields the nematic ordering field, a, equals 0.91 
( 5 CB-pSi 02 ) and 1.10 ( 6 CB-pSi 02 ) which in both cases cor¬ 
responds to a supercritical regime and appears in good agree¬ 
ment with recent optical birefringence studiesi^^. Whereas the 
changes of static capacitance characterises the behaviour of 
the effective (averaged) order parameter over the pore filling, 
analysis of the relaxation strengths of the slow and fast relax¬ 
ation processes provide an access to local molecular ordering 
in different parts of the pore filling. 

Thus, the molecular ordering in the surface or interface 
layer next to the pore walls and in the core region can be re¬ 
solved. The molecular ordering in the core region is reminis¬ 
cent of that in the bulk. It exhibits a strong increase in the ne¬ 
matic ordering just below the paranematic-to-nematic transi¬ 
tion point, TpN, and subsequent saturation at cooling. The sur¬ 
face ordering, on the other hand, evidently exhibits a contin¬ 
uous evolution with a smeared kink near Tp^ and asymptotic 
decreasing above this temperature. The molecular ordering 
in the core region of the pore filling influences the peripheral 
molecular layers, presumably via intermolecular interactions. 

Overall, the dielectric spectroscopy study on the calamatic 
nematic crystals 5CB and 6 CB confined to cylindrical nanochan¬ 
nels of monolithic silica membranes reveals features analo¬ 
gous to those found for 7CB in nanoporous silica substrates^. 
Therefore, we believe that the dielectric spectroscopy tech¬ 
nique along with the phenomenological approach and key phys¬ 
ical principles outlined here and in Ref.— are applicable to the 
important class of the cyanobiphenyl family and confined ne¬ 
matics in general. 

An inhomogeneous orientational behaviour has also been 
found in Molecular Dynamics simulations for rod-like LCs in¬ 
teracting with Gay-Berne potentials for confinement in nanochan- 
nel s^^i^^i^^ and in slit-pore geometryi^— . The qualitative re¬ 
sults of these studies agree with the findings presented here. 
Unfortunately, no partial fillings, in particular monolayer be¬ 
havior, have been explored in the simulation studies. This may 
be an interesting task for simulation studies in the future. 

Moreover, in agreement with experimental studie s 
molecular dynamics studies suggest a significant infiuence of 
interface roughness and/or micro structure on the anchoring 
strength and thus on the nematic-isotropic transition^^*^, an 
interplay which may be explorable in the near future given the 
availability of hierarchically tailorable porous solids^. 
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